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ABSTRACT 

We have carried out multi wavelength observations of the near-by {z = 0.046) rich, 
merging galaxy cluster Abell 3376 with the Murchison Widefield Array (MWA). As a 
part of the GaLactic and Extragalactic All-sky MWA survey (GLEAM), this cluster 
was observed at 88, 118, 154, 188 and 215 MHz. The known radio relics, towards the 
eastern and western peripheries of the cluster, were detected at all the frequencies. 
The relics, with a linear extent of ~ 1 Mpc each, are separated by ~ 2 Mpc. Gombining 
the current observations with those in the literature, we have obtained the spectra of 
these relics over the frequency range 80 ~ 1400 MHz. The spectra follow power laws, 
with a = —1.17 ± 0.06 and —1.37 ± 0.08 for the west and east relics, respectively 
{S oc z/“). Assuming the break frequency to be near the lower end of the spectrum 
we estimate the age of the relics to be ^ 0.4 Gyr. No diffuse radio emission from the 
central regions of the cluster (halo) was detected. The upper limit on the radio power 
of any possible halo that might be present in the cluster is a factor of 35 lower than 
that expected from the radio power and X-ray luminosity correlation for cluster halos. 
From this we conclude that the cluster halo is very extended (> 500 kpc) and/or most 
of the radio emission from the halo has decayed. The current limit on the halo radio 
power is a factor of ten lower than the existing upper limits with possible implications 
for models of halo formation. 

Key words: galaxies: clusters: individual (Abell 3376) - galaxies: clusters: intraclus¬ 
ter medium - radio halos - radio relics - MWA 


1 INTRODUCTION 

Galaxy clusters are some of the largest known gravitation¬ 
ally bound structures in the Universe composed primarily 
of dark matter along with galaxies and the X-ray emit¬ 
ting Intra-Cluster Medium (ICM) gas. Their masses are 
~ IO^^Mq while their sizes are ~ Mpc. According to the 
theory of hierarchical structure formation, galaxy clusters 
are formed by mergers of several smaller sub-clusters and 
groups driven by gravity. The process of cluster evolution 
through mergers is still ongoing and is one of the most ener¬ 
getic phenomena in the Universe with energies 10®® — 10®"* 


erg ( Hoeft et al.|2008 |. 

In some merging clusters Mpc-scale diffuse synchrotron 
radio emission, in the form of radio halos (in the central 
regions of the cluster) and/or radio relics (towards the pe¬ 
ripheries), have been detected. The radio halos have a typical 
surface brightness of a few mjy arcmin“® at 1.4 GHz (see 


Feretti et al. (20121 and Brunetti & Jones (20141 for recent 


reviews). Radio halos usually have a smooth morphology 
while relics are found to be elongated or arc-like in nature. 
The fraction of galaxy clusters which host such structures is 
rather small. Of all the known galaxy clusters below a red- 
shift of 0.2, ~ 5% contain radio halos and relics (Giovannini 
et aI.|[l9M |. This fraction increases to ~ 35% for higher X- 
ray luminosity clusters [Lx 10'*® erg s“*) (Venturi et al. 


20081 which have correspondingly higher radio power as well 
([Liang et al.||2000 Bacchi et al.||2003 jCassano et ar]|2006 


2007 


[Brunetti et al.||2007[ |2009| [Rudnick fc Lemmermanj 

2009). In a more recent study of galaxy clusters with red- 


shift 0.2 < < 0.4, Kale et al. (2013) estimated that 23% of 

galaxy clusters with X-ray luminosity Lx > 5 x 10'** erg s“* 
host radio halos. For clusters with Lx > 8 x 10** erg s“* 
the fraction increases to 31%. It has been observed, however, 
that highly disturbed, merging clusters are much more likely 
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to host halos and relics than relaxed clusters ( BuotepOOl 


Gassano et al. 2010). This led to the scenario that major 


cluster mergers could provide sufficient energy to the elec¬ 
trons in the ICM which, in the presence of existing magnetic 
field, produce synchrotron radiation that is detected as ha¬ 
los and relics. It should be mentioned here that even though 
both halos and relics are found predominantly in disturbed 
clusters, the processes responsible for their generation are 
thought to be different, as explained below. 

One of the currently accepted models for the genera¬ 
tion of radio halos is turbulent acceleration ( Petrosian |2001 


Brunetti et al.|2001 1 . During the process of a cluster merger. 


turbulence in the ICM gas re-accelerates the electron pop¬ 
ulation of the ICM to ultra-relativistic energies (7 > 1000, 
where 7 is the Lorentz factor). This turbulence also ampli¬ 
fies the magnetic field to a few pG (Garilli fc Taylor|[2002 


Subramanian et al. 2006). The synchrotron radiation from 


these electrons could be seen in the form of radio halos from 
the cluster centre. 

An alternative to the turbulent acceleration model is 


the hadronic (secondary) model ([Dennison[ 

[Blasi k\ 

Colafrancesco|1999 Pfrommer & Enfilin|2004 

1 , which posits 


that collisions between relativistic and thermal protons in 
the ICM are responsible for the generation of relativistic 
electrons. Proton-proton collisions result in the formation of 
pions. Neutral pions decay to produce gamma rays whereas 
charged pions decay to produce electrons and positrons. 
Since relativistic protons in galaxy clusters have a lifetime 
greater than that of the cluster itself, there will always be a 
fraction of the relativistic electrons in the ICM produced in 
such a manner. However, based on radio observations of mas¬ 


sive, X-ray luminous galaxy clusters, Brunetti et al. (2007) 


claimed that relativistic protons contribute less than a few 
percent of the energy in particles in galaxy clusters. Re¬ 
cently, observations of galaxy clusters at GeV energies were 


carried out using the Fermi and Cherenkov telescopes (Aha- 


[ronian et al.[[2009a|b| [Ackermann et aL[[ 2010 | [Aleksic et al. 
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2010). No gamma rays were detected and only upper lim¬ 


its to the emission were placed. Therefore, any contribution 
to synchrotron emission due to relativistic protons (indi¬ 
rectly) would be very small. It is believed that the level of 
radio emission produced in this process (secondary model) 
is about a factor of 10 below the radio emission produced in 
the turbulence model ( Brunetti fc Lazarian|20TT |. 

Radio relics, on the other hand, are created as a conse¬ 
quence of shock fronts travelling outwards from the cluster 
centre after a merger. The acceleration of electrons in the 
periphery of the cluster is best explained by the diffusive 
shock acceleration (DSA) theory ( [Blandford &; Eichlerjl987| 
Jones fc Ellison|1991 1. According to the theory, shock accel¬ 
eration is essentially a first-order Fermi acceleration process 
in which thermal particles are accelerated diffusively. They 
scatter elastically across the two fronts of a shock repeat¬ 
edly, gaining energy at each crossing of a front until they 
get accelerated to very high energies ( [Feretti et al.|[^12[ ). 
The best examples for the validity of this theory come in 
the form of a special class of relics, known as double relics, 
where two arc-like structures are seen on the opposite sides 
of the cluster centre. However there are very few of these 
double relics (e.g. Abell 3667 | Johnston- Hollitt|2003 1, Abell 
3376 ( Bagchi et al.|2006 ), CIZA J2242.8-1-5301 (van Weeren 


et al. 


2011| etc.) known. Rarer still are clusters that host 


both halos and double relics, even though both are thought 
to be a consequence of the same merger process (e.g. RXC 

CIZA J2242.8-b5301 


J1314.5-2515 (Feretti et al. 


2005) 


El 

Gordo (Lindner et al. 2014) etc.). The reason for the ab¬ 


sence of halos in merging systems is not well understood. 
One possibility is that the halo was created at the time of 
the merger but faded away over time. 

Radio halos and relics are systems where different phys¬ 
ical processes come into the picture. Detection or non¬ 
detection of diffuse radio emission depends on many fac¬ 
tors such as the mass ratio of the merging clusters, impact 
parameter of collision between the clusters, time since the 
merger, turbulence decay time scales, and synchrotron life¬ 
time of the radiating electrons ( Brunetti fc Jones|2014 ). In 
the case of radio relics, the inhomogeneity of the ICM also 
plays an important part. The dearth of detections of ra¬ 
dio halos and radio relics can be explained in many ways. 
Clusters with lower X-ray luminosity have correspondingly 
lower radio power and may not be bright enough to be de¬ 
tected by existing telescopes. Additionally, if one assumes 
the turbulent acceleration model that produces the radio 
halos, then the diffuse radio emission should display an ex¬ 
ponential cut-off at higher radio frequencies (> 1.4GHz). 
As a result most of the extended, diffuse emission will be 


fainter at higher frequencies (Cassano 2010). Furthermore, 


large angular extent diffuse emission is also likely to be re¬ 
solved out at higher frequencies. This is why it is necessary 
to observe and image galaxy clusters at low radio frequencies 
with greater surface brightness sensitivity than what GHz 
receivers can currently offer. Telescopes such as the MWA 
(Lonsdale et al. 2009[ Tingay et al. 2013| Bowman et al. 


2013), LOFAR (LOw Frequency ARray, van Haarlem et al. 


(2013)) and the upcoming SKA (Square Kilometre Array) 

I Dewdney et al. will be crucial in this endeavour. 

With the above motivations in mind, recently, |Hind-| 
[son et ah] ( |2014[ ) used commissioning data from the MWA 
to observe the galaxy cluster Abell 3667 at low frequencies 
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Figure 1. GMRT 325 MHz contours | |Kale et al.||20l'^ on the 
ROSAT X-ray (0.1-2.4 keV) broadband image. The units of the 
X-ray image are in counts s“^ arcmin“^. The contour levels start 
from 5.8 mjy beam“^ (Sfr) and increase by a factor of \/2. The 
resolution of the radio image is 39” X 39”. The X and the -|- 
represent the BCG2 and BCGl galaxies, respectively. The yellow 
solid line (drawn by the authors) marks the boundaries of the 
east and west relics. 


(120, 149, 180 and 226 MHz). A3667 has long been known 
to host two radio relics. Its northern most relic has been 


discussed in a number of papers spanning 45 years (Ekers 


1969[ Schilizzi fc McAdam|1975|[Goss et aI.|1982|[Rottgering 


et al. 1997), whilst the weaker second relic, though visible 


in the published 843 MHz images in Rottgering et al. (1997) 


was overlooked at that time, and not hrst mentioned in the 
literature until it was rediscovered in the 1.4 GHz data and 


the entire cluster re-imaged (Johnston-Hollitt et al. 2002 
Johnston-Hollitt]|2003 ) thus making A3667 the first double 


relic cluster. Imaging between 1.4 and 2.4 GHz has given 


rise to suggestions of a radio halo (Johnston-Hollitt et al. 


2002 Johnston-Hollitt 2003 Carretti et al. 2013), a possi¬ 


ble bridge connecting the two relics ( [Carretti et al.||2013[ ) 
or more recently a mini-halo off-set from the cluster centre 
( Riseley et al.|20f5 ). Assuming typical spectral characteris¬ 
tics, the putative halo and bridge were expected to be easily 
visible in the lower frequency MWA imaging. However, the 
MWA imaging did not confirm the detection of either struc¬ 
ture. Based on this MWA imaging it has been suggested 
that A3667 is seen in projection through a large diffuse spur 
of Galactic emission, which in turn might be responsible for 
these features ( [Hindson et al.||2014[ ). At the sensitivity of 
the MWA commissioning observations it was not possible to 
distinguish the Galactic features from a weak halo or bridge 
that might also be present, though strong limits were placed 
on the possible emission. In our current observations of the 
morphologically similar system, A3376, we expect to reach 
a sensitivity that is a factor of 5 better than the commis¬ 
sioning observations. 
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1.1 Abell 3376 


Abell 3376 is a rich, merging galaxy cluster at a redshift 


of 2 = 0.046 (Struble & Rood 19991. It has a n X-ray lu¬ 
minosity of I/xr o.i-2.4 kevi = 2.2 X 10"‘'*erg s“^ (Reipric h fc 


B6hririgeij2002 1 and a virial mass of ~ 3.64 x 10^*^ Mq (| Gi- 
rardi et al.|1998 1 . Optical observations | Escalera et al.|l994[ 
Durret et al.||20l3l indicate the merging status of the clus¬ 


ter by showing the two Brightest Cluster Galaxies (BCG) 
on the opposite ends of the cluster with the X-ray peak at 
BCG2 (MRC0600-399) (Fig. [^. The BCGl is a cD galaxy 
and is considered to be the centre of the main cluster. The 
BCG2 is an elliptical and its radio counterpart is a bent-jet 
galaxy. Machado & Lima Neto (20131 simulated the merger 


of A3376 and compared their simulation results with the X- 
ray morphology of the cluster. Based on these comparisons 
they suggest that the merger in A3376 occurred ~ 0.5 Gyr 
ago between subclusters having a mass ratio of 6:1. Very 


Large Array observations of this cluster at 1.4 GHz (Bagchi 


et al. 20061 showed the presence of two arc-like structures 
in the eastern and western peripheries of the cluster, ~ 1 
Mpc away from the cluster centre. These relics are assumed 
to be due to synchrotron emission from the ICM as a con¬ 
sequence of shock waves travelling outward from the cluster 
centre soon after the merger. Further radio observations at 
lower frequencies (150 and 325 MHz) using the Giant Meter- 
wave Radio Telescope ( Kale et al.||2012 l detected the relics 
but not a radio halo from the central regions of the cluster 
(Fig0. The east relic is found to be polarised up to 30% 
with aligned magnetic field vectors, whereas the west relic 
is polarised between 5 and 20% and the field vectors are not 
aligned in any particular direction ( Kale et al.||2012 l. 

The motivations for the current study are twofold. First, 
to detect the radio relics at low frequencies and establish 
their spectra over 80 - 1400 MHz. Second, to detect or put 
stringent upper limits on the radio power of a possible halo 
in the central regions of the cluster, which has eluded detec¬ 
tion so far. The paper is organised as follows: In section 2 
we discuss observations and data analysis. Section 3 contains 
the results and in section 4 we discuss the implications of 
these results. We conclude with the summary in Section 5. 
In this paper we adopt a A-GDM cosmology with Hq = 70 
km s“^ Mpc“^, Hm = 0.3 and Ha = 0.7. 


2 OBSERVATIONS AND ANALYSIS 

2.1 The GLEAM Survey 

The observations of A3376 were carried out as part of 
GaLactic and Fxtragalactic All-Sky MWA (GLEAM) Sur¬ 
vey (Wayth et ah, submitted). There are five observing 
frequency ranges: 72.3-103.04 MHz, 103.04-133.76 MHz, 
138.88-169.6 MHz, 169.6-200.32 MHz and 200.320-231.04 
MHz, each having a bandwidth of 30.72 MHz. The data 
were recorded with an integration time of 0.5 second and a 
spectral resolution of 40 kHz. The frequency range 134-138 
MHz was avoided due to the presence of satellite RFI in that 
range. 

The GLEAM survey was carried out in 
the drift-scan mode with the antennas pointing 
to one of the seven declination settings (<5 = 
-bl8.6°, -bl.6°, -13.0°, -26.7°, -40.0°, -55.0°, -72°) and 


letting the sky drift overhead. A single scan consists of 
observations for 112 seconds at each of the five frequencies 
available and then cycling through them. The data used for 
imaging A3376 in this paper were taken during November 
6-7, 2013 while covering the RA range 0-8h at declination 
S = —40°. This declination pointing is within a few 
arcminutes of the centre of the cluster (Fig. [^. 


2.2 Data reduction 


In principle, one could use the specialised pipeline software 
that has been set up to analyse the observations from the 
GLEAM survey (Hurley-Walker, in prep.). However, with 
the scope to try to analyse the MWA observations using 
a standard radio astronomy software we imaged the clus¬ 
ter at two frequencies viz. 154 and 215 MHz using GASA 
Q (Gommon Astronomy Software Applications). We chose 
these two frequencies because at 215 MHz the highest an¬ 
gular resolution in the GLEAM frequency range is achieved 
(151” X 145”) and at 154 MHz the MWA images obtained 
can be directly compared with the existing GMRT observa¬ 
tions at 150 MHz. The images at 154 and 215 MHz, pro¬ 
cessed with GASA, are further compared with the pipeline 
processed images at the same frequencies. For the remain¬ 
ing frequencies (88, 118 and 188 MHz) we use the pipeline 
processed images. 

The raw data from A3376 scans were converted to mea¬ 
surement sets readable by CASA using the COTTER soft¬ 


ware (Offringa et al. 20151. This software also performed 


some initial flagging and RFI removal using the AOFlagger 
library ( Offringa et al.|[2010 20121. The source 3C161 was 
used to create the normalised bandpass tables (task BAND¬ 
PASS in CASA) which were applied to all the A3376 scans. 
In all, a total of seven scans of A3376 were used - the tran¬ 
sit scan and three scans on either side of the source’s transit 
through the meridian - which covered a total time range of 
±30 minutes around transit at each frequency. The scans be¬ 
yond this time range were not considered due to the falling 
gain of the primary beam (i.e. the beam response was falling 
< 50%). Note that the full widths at half maxima of the 
primary beams of the MWA tiles are ~ 25° at 154 MHz 


I Hurley-Walker et al. 20141. The bandpass corrected mea¬ 


surement sets were then averaged in time and frequency to 
8 s and 400 kHz, respectively. Each scan was then individu¬ 
ally imaged and deconvolved in CASA using the task CLEAN. 
Each image produced (2048x2048 pixels) covered an area of 
34.13 X 34.13 deg^ in each axis with a pixel size of 1'. Decon¬ 
volution was performed using the Cotton-Schwab algorithm 
( Schwab! 19841 and widefield effects were taken into account 
using the w-projection algorithm ( Cornwell et al.||2008 l by 
setting the number of projection planes in CLEAN to 256. 

A variety of robust weighting schemes were tried to 
make images of A3376. In order to reach a balance be¬ 
tween sensitivity and resolution we used “Briggs” weight¬ 
ing I Briggs|1995 1 with robust value 0 for each scan at both 
polarisations. Note that Briggs weighting with values —5 
and ±5 correspond to purely uniform and purely natural 
weighting, respectively. The images were then sent through 
four rounds of phase self-calibration, at the end of which the 


^ http://casa.nrao.edu 
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phase errors had reduced to an rms of 0.2°. No further im¬ 
provement was seen with subsequent phase and amplitude 
self-calibration. The XX and YY polarisations were imaged 
separately for each scan. They were then corrected for pri¬ 
mary beam effects ( Sutinjo et al.|2014 l. 

Following the above procedure a total of seven images 
were produced for A3376 for each polarisation. These im¬ 
ages were examined separately and their rms values were 
estimated. Furthermore, the positions of unresolved sources 
in the GMRT 150 MHz image within a degree from A3376 
were compared with their respective positions in the seven 
MWA 154 MHz images. Across these seven images, the po¬ 
sitional offsets for all these sources were constant to within 
10% of the MWA synthesized beam size at 154 MHz. At 
each frequency, the individual primary beam corrected XX 
and YY images were weighted using the inverse variance of 
the rms in the respective images and then combined. The 
final XX and YY images were added together to produce 
the full total intensity (Stokes I) image. 

For the remaining frequencies (88, 118 and 188 MHz), 
we utilised the pipeline-processed GLEAM images, de¬ 
scribed in detail in Hurley-Walker et. al (in prep). For each 
scan, the basic steps of this pipeline are as follows, 

• Perform a first-pass calibration using a single bright 


source, applied to all observations (Hurley-Walker et al. 
[M4| ); 


Using WSGlean | Offringa et al. 20141, invert and clean 


multi-frequency-synthesis images in instrumental Stokes, 
across the full bandwidth, stopping at the first negative clean 
component; 


• Use the MWA primary beam model (Sutinjo et al. 20141 
to transform to astronomical Stokes; 

• Set Q, U and V to zero (i.e. assume that the sky is 
unpolarised) and use the same beam model to transform 
back to instrumental Stokes; 

• Use the sky model in a self-calibration loop to update 
the gain solutions; 

• Form new multi-frequency synthesis images. 

These individual images were then combined to produce 
the final total intensity images exactly as described earlier. 


2.3 Calibration 

The absolute flux density calibration of the MWA images 
was tied to those of the GMRT images at 150 and 325 MHz 
I Kale et al.|2012 1 . For GMRT calibration, 3G147 and 3G48 
as primary calibrators. The Perley-Butler (2010) scale was 
used for absolute flux density calibration. The secondary cal¬ 
ibrators used in the GMRT 150 and 325 MHz observations 
were 0521-207 and 0837-198, respectively. The flux densities 
of these two calibrators obtained from GMRT at the two fre¬ 
quencies agree to within 10% of their expected values based 
on the VLA calibrator list ( Kale et al.|2012 |. 

The MWA 154 MHz image was calibrated using the 
GMRT 150 MHz image. Unresolved sources in the GMRT 
150 MHz image within a degree from A3376 were used to 
estimate an average scaling factor between the GMRT 150 
MHz and the MWA 154 MHz images 

To calibrate the MWA 215 MHz image, the spectral in¬ 
dices of the same unresolved sources were estimated between 
the GMRT 150 MHz and 325 MHz images. The extrapolated 
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Figure 2. MWA 154 MHz (top) and 215 MHz (bottom) contour 
images. The 215 MHz contours are overlaid on ROSAT X-ray 
(0.1-2.4 keV) image. The units of the X-ray image are in counts 
s“^ arcmin“^. Both images have a robust weighting of 0. The 
154 MHz image has a resolution of 190” X 184" —85.6° while 
the 215 MHz image has a resolution of 151” X 145”, 68.2°. The 
synthesised beam is indicated at the bottom left corner in both 
the images. Contours for the 154 MHz and 215 MHz images start 
at 54 mjy beam“^ (Sc) and 33 mjy beam~^ (3fr), respectively, 
and increase by a factor of \/2. 


flux densities at 215 MHz were estimated for these sources. 
The MWA 215 MHz image was then scaled according to the 
average value of the scaling factor. The images at 88, 118, 
and 188 MHz were also calibrated in a manner similar to 
the one used for calibrating the MWA 215 MHz image. 

The images at 154 and 215 MHz made from GASA and 
from pipeline-processing were compared. At 154 MHz, the 
rms of the GASA processed image was 18 mJy beam“^ (Ta- 
ble[^ while that of the pipeline processed image was 13 mJy 
beam~^ . At 215 MHz, the rms of the GASA image was 11 
mJy beam“^ (Tablej^ while the pipeline image had an rms 
of 10 mJy beam“^ . The rms of the pipeline images could 
be expected to reduce further by mosaicking in narrower 
frequency channels, resulting in fewer excluded images and 
thus a higher signal-to-noise ratio (Hurley-Walker et ah, in 
prep). 
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Object Flux Density (mjy) 



88 MHz 

118 MHz 

154 MHz 

188 MHz 

215 MHz 

East Relic 
West Relic 
Central Region 

5476 ± 675 
2805 ± 348 
48.5 ±47 

3886 ± 577 
1926 ± 288 
59.5 ± 28.2 

2508 ± 465 
1351 ± 253 
22 ± 18 

1508 ± 169 
1172 ± 134 
20 ± 12 

2052 ± 350 
609 ±108 
17± 11 


Table 1. The integrated flux densities for the east and west relics along with their uncertainties are given in the first two rows. The 
third row gives the mean value of the flux density from the central region at different frequencies. These mean values represent upper 
limits to halo emission in the MWA synthesised beam (~160 kpc at 154 MHz). The errors quoted on these mean values for the central 
region are the rms values at the corresponding frequencies. 


East Relic 


West Relic 




Figure 3. Spectra for the east and west relics. These spectra include measurements from the GMRT observations at 150 and 325 MHz 
) |Kale et al.||201^ (blue), from the VLA observation at 1.4 GHz ( [Bagchi et al.|[2006^ (black) and from the current MWA observations 
at 88, 118, 154, 188, 215 MHz (red). Also displayed are the best fit lines for the respective data points as well as the spectral index 
estimates. 


3 RESULTS 

The MWA images of the Abell 3376 region at 154 and 215 
MHz are shown in Fig. The images at the other MWA 
frequencies are not presented here as they are essentially 
convolved versions of the above images. The Abell 3376 im¬ 
ages at 154 and 215 MHz show two prominent relics just 
outside the region of X-ray emission: one toward the east 
and the other toward the west of the X-ray emission (com¬ 
pare Fig.[^and[^. The two relics are separated by ~ 2 Mpc 
(~ 30') from each other. Both the east and west relics were 
detected at all the MWA frequencies. No diffuse radio emis¬ 
sion in the central regions of the cluster (halo) was detected 
down to the GLEAM sensitivity. Tablej^lists the integrated 
flux densities of the relics and the upper limits on the flux 
densities of a possible halo in the region between the relics. 

There are two primary sources of errors in the estima¬ 
tion of the flux densities of the relics in the MWA images. 
First, there is an error due to the uncertainties in the flux 
densities of the unresolved sources used for calibration of the 
MWA images. This error is estimated to be ~10%. Second, 
since the relics are extended sources in the MWA observa¬ 
tions, the errors in their flux density estimations will be the 
rms in the image multiplied by the square root of the ra¬ 
tio of the solid angle of the relic to that of the synthesized 
beam. Since these two sources of errors are unrelated, they 
are added in quadrature to estimate the final error on the 
flux densities of the relics that are quoted in Table 

The left and the right panels in Fig.j^show the spectra 
over the frequency range 80-1400 MHz, for the east and 


the west relics, respectively. The spectral index of the east 
relic is —1.37 ± 0.08 while that of the west relic is —1.17 ± 
0.06. A least-squares fitting method was used to estimate 
the spectral index values from the flux density values given 
in Table The reduced-x^ values for the fits on the east 
and west relics were 1.7 and 3.4, respectively. 


4 DISCUSSION 


Abell 3376 shows two arc-like radio relics on either side of 
the X-ray emission from the cluster. The arcs are oriented 
roughly perpendicularly to the merger axis of the cluster. 
Such a scenario is in agreement with the models that predict 
two outgoing shocks in the aftermath of the merger which 
are responsible for the radio emission seen from the arcs. 


According to the DSA theory (Blandford & Eichler 


19871 the slope of the power law spectrum of the number 


density of electrons, p, as a function of energy (N(E) oc E^), 
is related to the Mach number of the shock, A4, by 


+1 . 

P ^A42_i 


( 1 ) 


This includes the effects of particle ageing due to inverse 
Compton and synchrotron radiation losses under the as¬ 
sumption of continuous injection ( Sarazin|1999 Giacintucci 
|et al.|[^08[ ). The slope is related to the spectral index (a) 
of the synchrotron radiation by 
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P +1 


( 2 ) 


The spectral index is related to the flux density, S, by 
the standard relation, S' oc where v is frequency. For the 
east relic, the spectral index is estimated to be —1.37 ± 0.08 
which gives the Mach number 2.53±0.23. For the west relic, 
the spectral index is estimated to be —1.17 ± 0.06, which 
gives the Mach number 3.57 ± 0.58. 


Suzaku X-ray observations by Akamatsu et al. (20121 


estimate the Mach number of the East and West relics to be 
2.91±0.91. This number agrees with the value obtained from 


simulations (Machado & Lima Neto 20131 as well. The Mach 


numbers estimated here for both the relics are consistent 
with the X-ray estimates. 

The spectra of the relics show no breaks in the frequency 
range 80 to 1400 MHz. The shock acceleration is expected 
to produce an initial injection index for the synchrotorn ra¬ 
diation close to -0.8. Radiation losses induce a break such 
that the spectral index after the break steepens by 0.5 to 


-1.3 (See for example, the relic radio galaxies in Slee et al. 


(20011). The average spectral index of the two relics is 
-1.3. This is consistent with the break frequency being close 
to or lower than the lower end of the spectrum (~ 80 MHz) 
shown in Fig. The magnetic field in the cluster can be 
assumed to be ~ IpG ( Govoni et al.|20T3 |. We can estimate 
the spectral ages of the relics as, 


halos. The upper limit on the radio power of the A3376 halo 
is lower by a factor of ~35 compared to that expected from 
the best fit line. 

The upper limit on the radio power of the halo corre¬ 
sponds to a halo size of ~ 160 kpc (~ 3') at 154 MHz. This 
corresponds to the beam size at 154 MHz. If the halo is more 
extended than this, even with higher total radio power, it 
will be undetected because of limits on surface brightness 
sensitivity. For a canonical halo as extended as 500 kpc, the 
upper limit on the radio power would still be lower by a 
factor of ~ 4 as compared to that expected from the best-fit 
value. 

In order to explain why no radio halo is detected in 
this cluster we first need to discuss some of the time scales 
involved in the merger process. The separation between the 
two relics of A3376 is ~ 2 Mpc. Simulations by |Machad^ 


& Lima Neto (20131 estimate that the shock velocity, Ug, 
is ~ 2000 km s“^ with an estimated dynamical timescale, 
tdyn ~ 0.49 Gyr. This agrees with their simulations of the 
gas morphology which puts the cluster at approximately 0.5 
Gyr after central passage. 

MHD turbulence generated in cluster mergers has a de¬ 
cay time associated with it which is given by |Gassano &| 

poo^. 


Brunetti 


^-(Gyr)~ (^)~'. (4) 


t = 1060 




52 -b 52 


-[{l + z)u] 


(3) 


where B is the magnetic field in the source, Bjc is the effec¬ 
tive Inverse Gompton magnetic field, 2 is the redshift of the 
source and v is the break frequency ( Slee et al.|2001 1 . Based 
on these values, the spectral ages of the relics is ~ 0.4 Gyr. 

We have made the assumption here that the break fre¬ 
quency is 80 MHz. In reality, this break could be lower than 
this which would imply that the spectral age is greater than 
0.4 Gyr. This estimated spectral age of the relics is, never¬ 
theless, consistent with the age of the cluster (~ 0.5 Gyr) as 
obtained from simulations discussed later in this section. 

The cluster merger responsible for the generation of the 
observed relics also generates turbulence in the ICM which 
is believed to be able to re-accelerate electrons in the region 


between the relics (Petrosian 2001 Brunetti et al. 20011. 


The region between the relics is then expected to host a radio 
halo. No such radio halo is detected at the MWA frequencies. 
However, we have been able to put an upper limit to the flux 
density from a possible halo in the central region of A3376. 
Table [ 1 ] lists the mean flux densities observed in the central 
regions of A3376 in the MWA images. The mean value of the 
spectral index of halos is —1.34±0.28 (Feretti et al. ( 2012[ )). 
An upper limit to the radio power of a halo at 1.4 GHz can 
be estimated using the upper limit at 154 MHz (Table 
and extrapolating it to 1400 MHz using the mean value of 
the spectral index quoted here. This power corresponds to 
6 X 10^^ W Hz-b 

The upper limit is plotted in Fig where the radio 
powers at 1.4 GHz and the corresponding X-ray luminosities 
of known halos and relics are also plotted ( [Brunetti et al.| 
2009). The solid line represents the best fit to the known 


where Ui is the relative velocity of the merging clusters at the 
time of impact and Tjt is defined as the fraction of energy in 
the form of turbulence that is converted into magnetosonic 
(MS) waves. This means that the value of rjt defines the effi¬ 
ciency with which the electron population is re-accelerated. 
The value of rji can be approximated to 0.25 by using the 
constraint that the slope of the radio spectrum of the halos 


is in the range, a = 1.1 — 1.5 (Kempner & Sarazin 2001). 


The value of Vj i s estimated from simulations ([Machado & 
Lima Neto '^01 oh +„ u„ .. icon 1 ,™ „-i 


20131 to be ~ 1500 km s . The typical injec¬ 
tion length scale (Linj) of turbulence in galaxy clusters is 


200 — 300 kpc ( Brunetti fc Lazarian[[20TT ). With these val 
ues in the above equation, the decay time of turbulence for a 
cluster like A3376 is estimated to be 0.25 — 0.4 Gyr. Consid¬ 
ering that it is now ~ 0.5 Gyr since core passage of the clus¬ 
ter, most of the turbulence that would have been generated 
at the time of merger would have dissipated. Furthermore, 
the half-life of these radio emitting electrons is short, ~ 0.3 
Gyr ( [Slee et S][2001[ ). Any radio halo that was generated 
during the cluster merger is likely dissipated. 

It has been postulated that in galaxy clusters there 
should be a long-lived component of cluster-wide diffuse ra¬ 
dio emission that is generated by relativistic electrons which 
are continuously injected into the IGM by relativistic pro¬ 
ton - thermal proton collisions ( Brunetti fc Lazarian|20lT |. 
The radio power of diffuse emission expected from these sec¬ 
ondary particles is expected to be at a level that is about an 
order of magnitude lower than that of the radio powers in 
halos that are detected. The upper limit on the radio power 
of a halo in Abell 3376 is nearing such levels. While the cur¬ 
rent limit is close to such a level, further lowering of these 
limits either in this cluster or in other clusters will have 
implications for the secondary models of halo production. 
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Figure 4. Plot showing the Lx — Pi .4 correlation for halos and relics. The filled circles represent all known radio halos. The empty circles 
represent radio relics jFeretti et al.|2012[|. The solid line represents the best fit for the Lx — Pi .4 relation for halos (log(Pi, 4 ) — 24.5 = 
0.195 + 2.06[log(Ljy) — 45], [Brunetti et al.|p009[ l). The arrows represent upper limits on halo emission ( [Venturi et al.|2008| [Kale et al.] 
|2013| |. The cross represents the A3376 relics. The diamond symbol shows the upper limit to halo emission in A3376 as estimated by 
GMRT while the arrow below it represents the new upper limit as estimated by the MWA. 


5 SUMMARY 

In this paper, we present low frequency radio observations 
of Abell 3376 using the Murchison Widefield Array at 88, 
118, 154, 188 and 215 MHz. The observations were carried 
out as part of the GaLactic and Extragalactic All-sky MWA 
(GLEAM) survey. The observations were analysed using 
Gommon Astronomy Software Applications (CASA) tools 
as well as specialised software packages written to process 
the GLEAM data. Images were produced at all frequencies. 
Towards the east and west peripheries of the cluster bright 
arcs of radio emission (relics) were detected at all frequen¬ 
cies. These arcs, with a linear extent of ~ 1 Mpc each, are a 
result of outgoing shocks produced by the merger event and 
are separated by ~ 1 Mpc from the cluster centre. Spectra 
for both the relics were estimated using the flux densities es¬ 
timated from the MWA images along with those from previ¬ 
ous observations of GMRT at 150 and 325 MHz and of VLA 
at 1400 MHz. The spectral indices of the east and west relics 
over the frequency range 80-1400 MHz are —1.37 ± 0.08 and 
— 1.17 ± 0.06 respectively. The Mach numbers of the shocks 
estimated from these indices are consistent with those esti¬ 
mated from Suzaku X-ray observations. Assuming that the 
break frequency due to radiation losses corresponds to the 
lower end of the spectrum, the age of the relics can be esti¬ 
mated to be ~ 0.37 Gyr. This age is consistent with the age 
of the cluster estimated from simulations to be ~ 0.5 Gyr. 

No diffuse radio emission was detected in the region be¬ 
tween the relics. The upper limit on the radio power of any 
possible halo of 160 kpc size is ~ 6x 10^^ W Hz“^ at 1.4 
GHz. This is a factor of 35 lower than that expected from 
the correlation between radio powers and X-ray luminosi¬ 


ties of cluster halos. It is very likely that the radio emission 
from the halo has substantially decayed over its age (~ 0.37 
Gyr) considering that this time scale is comparable to the 
life time of the radiating electrons at the MWA observing 
frequencies. The upper limit on the radio power of a halo 
in Abell 3376 is a factor of 10 lower than most upper limits 
available so far and is nearing levels at which radio emission 
from secondary electrons produced in relativistic proton - 
thermal proton collisions is expected. Given the uncertain¬ 
ties in the hadronic model | Brunetti fc Lazarian|20TT I it is 
not possible to discuss the validity of the model based on 
the observations of the cluster alone. However, observations 
of a large number of clusters with MWA and LOFAR at the 
current limits of detection or using the SKA at improved 
sensitivities will have implications on the hadronic model. 
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